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ABSTRACT 

We investigate properties of the correlation function of clusters of galaxies using 
geometrical models. We show that the correlation function contains useful information 
on the geometry of the distribution of clusters. On small scales the correlation func- 
tion depends on the shape and the size of superclusters. On large scales it describes 
the geometry of the distribution of superclusters. If superclusters are distributed ran- 
domly then the correlation function on large scales is featureless. If superclusters have 
a quasi-regular distribution then this regularity can be detected and measured by the 
correlation function. Superclusters of galaxies separated by large voids produce a cor- 
relation function with a minimum which corresponds to the mean separation between 
centres of superclusters and voids, followed by a secondary maximum corresponding 
to the distance between superclusters across voids. If superclusters and voids have a 
tendency to form a regular lattice then the correlation function on large scales has 
quasi-regularly spaced maxima and minima of decaying amplitude; i.e., it is oscillating. 
The period of oscillations is equal to the step size of the grid of the lattice. 

We also calculate the power spectrum and the void diameter distribution for our 
models and compare the geometrical information of the correlation function with other 
statistics. We find that geometric properties (the regularity of the distribution of 
clusters on large scales) are better quantified by the correlation function. 

We also analyse errors in the correlation function and the power spectrum by 
generating random realizations of models and finding the scatter of these realizations. 

Key words: cosmology; clustering - large-scale structure of the universe - methods; 
numerical 



1 INTRODUCTION 

The correlation function is one of the most frequently used 
statistics in cosmology. Conventionally the correlation func- 
tion and the power spectrum are considered together to char- 
acterise the basic clustering properties of galaxies or of clus- 
ters of galaxies. They complement each other in the sense 
that the power spectrum describes the structure in Fourier 
space and the correlation function in real space. 

Intuitively, it is clear that, on small scales compara- 
ble to the size of systems of galaxies, the galaxy correlation 
function should characterise the distribution of galaxies in 
systems of galaxies. On these scales the galaxy correlation 
function is large - this is the manifestation of the clustering 



of galaxies. With increasing scale the clustering gets weaker 
and at some scale the correlation function approaches zero. 
Traditionally the correlation function is expressed in log- 
log units and is calculated only up to the zero crossing sep- 
aration. This separation is sometimes used as a parameter 
characterising the transition from clustering to homogeneity. 

The correlation function of clusters of galaxies has a 
similar behaviour, only it is shifted towards larger scales. 
The clustering scale is here about a factor of 4 times larger 
than that for galaxies. 

On large scales the correlation function of galaxies de- 
pends primarily on the distribution of systems of galaxies. 
Thus, properties of the distribution of these systems influ- 
ence the behaviour of the correlation function on large scales. 
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Similarly, the correlation function of clusters of galaxies on 
large scales depends on the distribution of systems of clus- 
ters, i.e. superclusters of galaxies. Einasto et al. (1997a,b, 
hereafter E97 and Paper II) determined the correlation func- 
tion of clusters belonging to rich superclusters using data in 
a distance interval ~ 650 h^ 1 Mpc. The main finding of this 
study was the detection of a series of almost regularly spaced 
maxima and minima. This property of the correlation func- 
tion was called oscillations. The period of oscillations was 
found to be 115 ± 15 h' 1 Mpc. 

This behaviour of the correlation function emphasises 
the presence of certain regularities in the distribution of 
clusters of galaxies. A regular network of clusters of ap- 
proximately the same scale was actually noticed by Tully 
et al. (1992) and Einasto et al. (1994, 1997c, hereafter Pa- 
per I). 

These considerations have motivated us to perform an 
analysis of the correlation function as a geometrical statistic. 
Such analysis helps us to understand the relation between 
the geometry of the distribution of clusters and superclus- 
ters and the correlation function. We shall use simple ge- 
ometrical models of the particle distribution and calculate 
the correlation function for these models with the aim to 
find the relation between the correlation function and the 
geometry of the cluster distribution. Geometric information 
contained in the correlation function has been studied so far 
only for relatively small separations. Bahcall, Henriksen, & 
Smith (1989) used a model in which galaxies were located 
on the surfaces of spherical shells and clusters of galaxies 
resided at the intersections of the shells. Several models are 
discussed by Einasto (1992) with the aim of giving a geo- 
metrical interpretation for the galaxy correlation function on 
small scales. It was also noted that if the systems of galaxies 
are located in a regular grid then the correlation function 
has a secondary maximum at a scale which corresponds to 
the grid length (step). 

In contrast to previous studies we focus our attention 
in the present paper on large scales. On these scales the er- 
rors of the correlation function are of the same order as the 
function itself, thus a careful error analysis is crucial. Ad- 
ditionally we found for our models the power spectrum and 
the distribution of void radii in order to compare the rela- 
tionship of these statistics to sample distribution geometry. 

The paper is organised as follows. We start the paper 
with a short description of the basic formulae and meth- 
ods we use in calculations (Section 2). Then we use analytic 
models for the spectra and calculate respective correlation 
functions to get a general picture of the behaviour of the 
correlation function on large scales (Section 3). Thereafter 
we describe geometrical models of the distribution of clus- 
ters used to investigate properties of the correlation func- 
tion (Section 4). The discussion of the correlation function 
for these models, the comparison of this function with the 
power spectrum and void radius statistics, and the discus- 
sion of the influence of the sample shape follow (Section 5). 
Next we discuss errors of the correlation function (Section 
6). Finally we discuss geometric properties of the correla- 
tion function on small scales (Section 7). The paper ends 
with the main conclusions. 



2 BASIC FORMULAE 

The two-point correlation function £ (r) is defined as the ex- 
cess over Poisson of the joint probability of finding objects 
in two volume elements separated by r and averaged over 
a very large volume (Peebles 1980). It is possible, in prac- 
tice, to determine this function only for a limited volume, 
which is called the estimate of the correlation function. We 
shall use the term "correlation function" for its estimate, 
and calculate it using the classical formula: 



Cl 1 (RR(r)) n 2 



(1) 



where (DD(r)} is the number of pairs of galaxies (or clusters 
of galaxies) in the range of distances r ± dr/2, dr is the bin 
size, (RR(r)) is the respective number of pairs in a Poisson 
sample of points, n and ur are the mean number densities 
of clusters in respective samples, and brackets (. . .) mean 
ensemble average. The summation is over the whole volume 
under study, and it is assumed that the galaxy and Poisson 
samples have identical shape, volume and selection function. 

We define our model samples in cubic volumes of side- 
length L, and calculate the correlation function for the whole 
distance interval from r = to r = L, dividing this interval 
linearly into 100 equal bins. For a given sample shape and 
selection function we calculate (RR(r)) /n 2 R as a mean value 
of ten realizations of random samples of 2000 clusters with 
the same selection function as the model sample. In this 
paper we use only simple selection functions: in addition 
to cubical samples we consider double-conical samples (see 
below). 

To suppress random noise we smooth the correlation 
function with a Gaussian window. After smoothing the Pois- 
son error is virtually absent and the residual error is deter- 
mined by the cosmic variance of the geometry of the particle 
distribution, i.e. by the variation of the correlation function 
found for different volumes of space. 

We calculate also the power spectrum, P(k), for our 
models. For this purpose we can use the relation between 
the power spectrum and the correlation function, as they 
form a pair of Fourier transformations: 
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Here the wavenumber k is measured in units of h Mpc -1 , and 
is related to the wavelength, A = 2n/k. To find the power 
spectrum from the correlation function we could use the 
Eq. (3). However, it is more convenient to use direct meth- 
ods which are similar to the calculation of the spectrum of 
./V-body models of structure evolution, calculating first the 
distribution of clusters in the Fourier space and then aver- 
aging over various wavenumber directions. The distribution 
of clusters in fc-space is found with the fast Fourier trans- 
form using 64 3 cells, which allows us to find 32 modes in 
the spectrum. This method assumes that the whole space is 
filled with similar periodic cubes. 

In addition to geometric models where the distribution 
of clusters in real space is given initially, we also find for com- 
parison the relation between the correlation function and 
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Figure 1. The power spectra (left) and correlation functions (right) for a double power model with sharp transition (dotted line) and 
with smooth transition (dashed line) and for the standard CDM model (solid line). 



power spectrum analytically for some simple power spec- 
trum models. In this case we use the Eq. (2). 



3 ANALYTIC MODELS 

Following Eq. (2) we have calculated the correlation function 
for several power spectra. In general, for any given power 
spectrum the correlation function cannot be calculated an- 
alytically. However, we can determine £(r) using a simple 
ansatz for P(k): 



P(k) 



Ak n , k < k ; 

Ak%(k/ko) m , k > k , 



(4) 



where n = 1 is the power index on large wavelengths 
(Harrison-Zeldovich spectrum), the negative m is the index 
on small wavelengths, fco is the transition wavenumber, and 
A is the normalisation amplitude. From observations m is 
expected to be about —1.5. Thus, we have to integrate 
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where x = kr, and y — kor depends on the place fco of 
the maximum of the power spectrum. The solution of the 
integral yields 
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where S denotes Fresnel's sine-integral. Obviously, the cor- 
relation function £(r) is oscillating. The amplitude scales 
with fco and the radial dependence with kor. The integrals 
in Eq.(5) (the power spectrum left and right of the max- 
imum) lead to approximately the same contribution in £, 
however with opposite signs. The resulting oscillations are 
rapidly damped (oc r -3 ). 

As a second model we consider a smooth transition 
between positive and negative spectral indices (Peacock & 



West 1992, Einasto & Gramann 1993) 

Ak n 
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The third model is the standard cold dark matter model 
(SCDM) with the transfer function of Bond and Efstathiou 
(1984), where F = flh — 0.5. The transition between positive 
and negative spectral indices is slower than in the second 
model. The maximum of the first two models is chosen to 
be at fco = 0.0562/i~ 1 Mpc in accordance with the SCDM 
model. The amplitude is taken in accordance with COBE 
data by adopting A = 6.4 x 10 B /i~ 4 Mpc 4 (White, Scott & 
Silk 1994). 

A rapid and precise method to calculate numerically the 
integral (2) is given by Press et al. (1992). The dispersion of 
density fluctuations is equal to the value of the correlation 
function at zero separation, 



p(k)k 2 dk. 



(8) 



In all cases of practical interest the spectral index is m > —3; 
i.e. the integral (8) is not convergent. Since the dispersion of 
density fluctuations has a finite value, we conclude that on 
small scales the real spectral index must decrease so that the 
spectrum becomes steeper. The physical reason for this be- 
haviour is galaxy formation; i.e. the power spectrum arises 
from objects of finite size. This effect can be formally de- 
scribed by an exponential cutoff of the spectrum with an 
effective wavenumber k s — 50 h Mpc -1 corresponding to 
a scale w 0.1 h^ 1 Mpc. This parameter does not change 
the correlation function on scales of interest for the present 
study. The spectra and respective correlation functions are 
shown in Figure 1. 

A smooth transition between positive and negative 
spectral indices can be approximated by a number of sec- 
tions with constant spectral index m . The index changes at 
wave-numbers fc;. Then the correlation function is a super- 
position of given analytic periodic functions with different 
periods 2n/ki, so the resulting correlation function does not 
oscillate (Gottlober 1996). 

It is evident that the amplitude of the correlation func- 
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Figure 2. The power spectra (left panels) and correlation functions (right panels) for the random superclustcr model (upper panels) and 
the Voronoi model (lower panels). Solid curves indicate the mean function for ten realizations, dashed lines show the la error corridor 
calculated from the scatter of individual functions. 



tion is proportional to the amplitude of the power spec- 
trum. Furthermore, in the first example with an analytical 
solution, the scaling of the correlation function with ko is 
obvious. In general the power spectrum can be written as 

P(k) = AkT 2 (<;k), (9) 

where A is the normalisation constant, £ defines the place 
of the maximum (£ = k^ 1 in the ansatz (4) and £ « 
20(f2/! 2 )" 1 Mpc for the CDM model, and T is the trans- 
fer function. Thus the correlation function is a function of 

x = r/C, 

£{x) = C*f{x), (10) 

where 

fix) = - [ T 2 (y) S m(yx)dy. (11) 
x Jo 

All correlation functions scale in the same manner as the 
analytic solution (6). 

Let us finally consider a correction Pi(k) to the power 
spectrum P(k) = Po(k) + Pi(k) which leads to £(r) = 
(,o(f) + £i(f)- Obviously, this correction leads to superim- 
posed oscillations of the correlation function if it has the 
form of the power spectrum Eq. (4) (see Section 4.3). The 
addition of such a correction with a relative power of about 
25 % leads already to significant oscillations in the correla- 
tion function. 



4 GEOMETRICAL MODELS OF THE 
LARGE-SCALE STRUCTURE 

In this Section we describe and analyse the correlation func- 
tions and the power spectra for a number of geometrical toy 
models of the distribution of clusters of galaxies. Models are 
geometrical in nature; i.e. in these models not the power 
spectrum is given initially but the distribution of clusters. 
Models can be divided into three basic types: random mod- 
els, regular models, and mixed models. In all models we ex- 
plicitly generate clusters that in most models are clustered 
into superclusters, and designate the models correspond- 
ingly by capital letters CL or SC. The following small letters 
in the model designation indicate the character of distribu- 
tion; there follow two digits which indicate the dispersion 
of the location of superclusters around the mean position. 
These designations shall be described in detail below. 

For each model we generate ten random realisations, 
calculate the mean smoothed correlation function, and esti- 
mate the parameters of the correlation function as described 
below. We find also the local rms deviations of individual 
correlation functions from the mean value. These deviations 
define the error corridor for the correlation function. 



4.1 Random models 

Random models are designated by "ran" in the Tables and 
Figures. We have generated random cluster, supercluster, 
and void models. 
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Figure 3. The power spectra and correlation functions for the regular net (upper panels) and regular rod models (lower panels). 



Random cluster samples are generated with the goal 
of checking the formulae for the calculation of errors. The 
number of clusters was taken to be from several hundred to 
several thousand. 

Random superclusters. In this model supercluster cen- 
tres are positioned randomly in rectangular coordinates. The 
multiplicity of the supercluster (the number of clusters in 
supercluster N ct ) was chosen randomly in agreement with 
the multiplicity distribution of real superclusters, as derived 
in Paper I. The multiplicity distribution of real superclus- 
ters is from 1 (isolated clusters like the Virgo cluster) to 32 
(the Shapley supercluster). For simplicity only superclusters 
of richness N c i = 1,2,4,8,16, and 32 are generated. Mem- 
ber clusters of superclusters of richness 2 and higher are lo- 
cated around the centre of the supercluster randomly with 
an isothermal radial distribution. Actual superclusters are 
elongated and we have checked the influence of the elonga- 
tion on the correlation function. In a model with elongated 
superclusters negative sections of the correlation function 
have smaller amplitudes. Positions of maxima and minima 
are not influenced. The radius of superclusters was chosen as 
a slow function of the supercluster richness, Rsc = RoN^ 25 . 
The parameter Ro determines the cluster correlation radius; 
we used a value 15 h^ 1 Mpc. These parameters determine 
the cluster correlation function and the power spectrum on 
small scales; values chosen yield good approximations to ob- 
served functions. In a cubic volume of size 700 h' 1 Mpc 1200 
superclusters are generated, the total number of clusters be- 
ing a random number around 3000. In a sample of randomly 
located superclusters some .superclusters are located closely 
together and merge. To find the actual multiplicity distri- 



bution of superclusters a clustering algorithm was applied, 
as for real cluster samples using a neighbourhood radius 24 
h~ x Mpc. For details of supercluster definition see Paper I. 

Random voids of diameter 100 hT 1 Mpc are distributed 
randomly and clusters of galaxies are located randomly out- 
side voids, i.e. no superclusters are generated. This diameter 
is equal to the mean diameter of voids in samples of rich 
clusters of galaxies (Einasto et al. 1994, Paper I). 

The Voronoi tessellation model is a variant of the ran- 
dom void model. In this case centres of voids have also been 
located randomly, but in contrast to the previous model 
clusters are not distributed randomly outside voids. Instead, 
void centres are considered as seeds, and clusters of galaxies 
are located in corners of a structure which is formed by ex- 
panding volumes of seeds. Superclusters in this model form 
by clustering of clusters if their mutual distances are smaller 
than the neighbourhood radius 24 h~ x Mpc used in the def- 
inition of superclusters. The Voronoi model has been inves- 
tigated by a number of authors (van de Weygaert & Icke 
1989; Ikeuchi & Turner 1991; Williams, Peacock & Heavens 
1991; SubbaRao & Szalay 1992; van de Weygaert 1994). We 
designate the model by "vor" . It has one free parameter, the 
number density of void centres. This determines the mean 
separation of superclusters across voids which is taken as the 
characteristic length of the model. We generated ten models 
of box size 700 h^ 1 Mpc with 431 void centres, in which case 
the characteristic scale is 115 hT 1 Mpc, a value observed in 
our study of superclusters and voids in Paper I, and in the 
study of the cluster correlation function in Paper II. 

Results for the correlation functions and power spectra 
of random models are shown in Figure 2. All random mod- 
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els have a correlation function with no sign of oscillations. 
The correlation function of the random supercluster model 
approaches zero on large scales immediately after the initial 
maximum. The correlation function of the Voronoi model is 
different. It has a minimum at a separation between ~ 60 
and ~ 100 hT x Mpc and one secondary maximum at ~ 140 
h^ 1 Mpc. We discuss differences between the random super- 
cluster and Voronoi models later. 

The analysis of the correlation function of real cluster 
samples has shown that it depends strongly on the richness 
of superclusters (see Paper II). We have thus made separate 
analyses of the correlation function for our random models 
for clusters located in rich and poor superclusters. We shall 
discuss this problem in more detail below (Section 7). 

4.2 Regular models 

In these models clusters or superclusters are located, with 
a certain degree of fuzziness, on the corners or on the 
edges between corners (rods) of a regular rectangular three- 
dimensional grid. The fuzziness of the model was realized 
by adding to positions of clusters or superclusters random 
shifts around the exact position on the corner or the edge 
of the grid. Depending on the degree of the fuzziness we 
can speak of strictly regular models with no fuzziness or 
of quasi-regular models with variable fuzziness (see also 
Einasto 1992). 

Clusters on regular rods. In this model clusters of galax- 
ies are located randomly along rods, regularly spaced in a 
rectangular grid with a step 115 h^ 1 Mpc. Rods are ran- 
domly positioned around the exact values defined by the 
regular grid with a scatter ±20 h~ Mpc, which is equal to 
the scatter of mean diameters of voids defined by rich clus- 
ters of galaxies (Einasto et al. 1994, Paper I). In our nam- 
ing convention, the model is designated by "rod". Rods are 
cylindrical, and the diameter of cylinders was taken equal to 
be the mean minor diameter of superclusters, ~ 20 h^ 1 Mpc 
(Einasto et al. 1994, Jaaniste et al. 1997). Rods cross at the 
corners of the grid, and at these points density enhancements 
form which can be considered as superclusters. 

Superclusters on regular rods. In this model instead of 
clusters, superclusters are located randomly along regularly 
spaced rods with step 115 ± 20 h~ x Mpc. 

Superclusters in corners of a regular grid. In this case 
superclusters are located only in corners of a cubic lattice 
with constant grid step. In our naming convention, this 
model is designated as "cor" . Clusters of galaxies in super- 
clusters are generated as in the random supercluster model. 
Random shifts to supercluster positions around corners are 
added with different values of the mean shift. 

Superclusters in a regular net. This model is a superpo- 
sition of the model with superclusters located in corners of 
the net, and additional random clusters located along regu- 
lar rods; this model is designated as "net". Clusters in rods 
represent isolated clusters and cluster filaments which are 
concentrated mostly between superclusters. Here, too, vari- 
able random shifts to mean positions in the regular net can 
be added. 

The correlation functions and power spectra for regular 
models are shown in Figure 3. 

The extreme case of regularity is a model with super- 
clusters of equal multiplicity which are positioned exactly in 



corners of the regular grid. The only free parameter it has 
is the step size of the grid. The correlation function and the 
spectrum consist of a number of isolated peaks correspond- 
ing to mutual distances of all possible pairs of distances of 
corners of the grid. The first peak of the correlation function 
is located in a separation equal to grid step, the next one 
to the shortest diagonal, and so on. There exists no regular 
oscillations of the correlation function of the type observed 
in analytical models with sharp maximum in the spectrum. 
This example shows that a strictly regular structure cannot 
produce regular oscillations in the correlation function. 

If we add some flexibility to the regular corner model 
by adding random shifts to supercluster positions then 
neighbouring peaks of the correlation function join to form 
broader maxima. Adding clusters on rods makes the cor- 
relation function even more smooth. One such example is 
shown in upper panels of Figure 3. If we make the geom- 
etry of the model even more random and allow a random 
position of superclusters on rods we get a correlation with 
even more regular oscillations, as seen in the lower panels of 
Figure 3. This example shows that in order to get a corre- 
lation function with regular oscillations there must be a lot 
of freedom in the supercluster positions. These models are 
regular only in a restricted sense (approximately equal step 
between rods). 

With changing random shifts the amplitude of oscilla- 
tions also changes. It is larger the more accurately super- 
clusters are positioned in corners or rods of the regular net. 
The scatter of superclusters around the mean position used 
in our toy models is approximately equal to the scatter of di- 
ameters of voids between real superclusters. Thus our mod- 
els should have correlation functions with an amplitude of 
oscillations which is close to the amplitude of oscillations 
of a population of real clusters, if clusters form a regular 
network of superclusters and voids with this scatter. The 
analysis made in Paper II shows that this is indeed the case. 

4.3 Mixed model 

The analysis of the correlation function of real cluster sam- 
ples has shown that clusters in high- and low-density envi- 
ronment have rather different correlation functions. Clusters 
in high-density environment (the members of very rich su- 
perclusters) have an oscillating correlation function whereas 
clusters in low-density regions (members of poor superclus- 
ters and isolated clusters) have an almost zero correlation at 
large separations (Paper II). To model this behaviour of the 
correlation function we have constructed a mixed model. 

The model has two populations of clusters: rich super- 
clusters are located randomly in regular rods and form the 
quasi-regular population, poor superclusters are located ran- 
domly throughout the whole volume of the model. For sim- 
plicity we use here randomly located superclusters as a low- 
density cluster population, taking into account the fact that 
the correlation function is not sensitive to details of the dis- 
tribution of clusters in low density environments (real clus- 
ters of the low density population lie near void walls, and 
void interiors are empty of rich clusters, see Paper I). One 
third of all clusters belong to the quasi-regular population, 
two thirds to the random one. 

The correlation function and spectrum for this model 
is presented in Figure 4. The correlation function is oscil- 
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Figure 4. The power spectra (left panels) and correlation functions (right panels) for the mixed model SC. mix. 15 and the double 
power-law model DPS. 6, shown in the upper and lower panels, respectively. 



lating but with a very low amplitude. The whole observed 
cluster sample has a similar correlation function, and the 
amplitudes of oscillations are also similar (Paper II). This 
example shows that our mixed model is a good approxima- 
tion to the real cluster sample. 

The mixed model is compared with a iV-body model of 
structure evolution based on the double power-law spectrum 
with a sharp transition between short and long wavelength 
spectral indexes, Eq. (4). The model was calculated similarly 
to the CDM model and was also used to estimate the cos- 
mic variance of the correlation function (see Section 6). The 
model is designated as DPS. 6. The oscillatory behaviour of 
the double power-law model is less regular; the first mini- 
mum of the correlation function has much larger amplitude. 
Amplitudes of maxima are approximately as high as in the 
case of the mixed model, the mean period of oscillations is 
also similar, but the scatter of individual periods is larger. 



4.4 Parameters of oscillations of the correlation 
function 

We shall use the following parameters to characterise oscilla- 
tions of the correlation function: the mean separation of the 
first minimum from zero, r m i„; the position and amplitude 
of the first secondary maximum, r max and A max ; differences 
between the subsequent maxima, Ay = r maxi — r max j\ and 
the mean value of differences, derived from the first 4 dif- 
ferences between maxima and from the first 4 differences 
between minima, A mean . Scaling parameters can be used to 
derive the value of the true period of oscillations, P, which 



we take to be equal to the grid step of the net of the quasi- 
regular model, and to the wavelength of the maximum of 
the spectrum, Ao = 27r/fco. 

Values of these parameters for models with at least one 
secondary minimum and maximum are given in Table 1. 
They are derived from the mean smoothed correlation func- 
tion of the particular model, calculated from ten realizations. 
To avoid the decrease by smoothing the amplitude has been 
estimated from unsmoothed data. 

This Table shows that these parameters vary only 
within rather narrow limits. The position of the first sec- 
ondary maximum is almost the same for all geometrical 
models with identical grid step (and period P). Differences 
between maxima and minima are larger, but not much. 
These differences are rather systematic: in all our geometric 
toy models A21 > P, and A32 < P- The characteristic scale 
of the supercluster-void network, P, can be calculated from 
the above parameters using the mean values of the following 
parameters: 



fi = 



1.20 ±0.05; 



and 



f 2 i = = 1.06 ±0.05; 



f 3 2 = % = 0.93 ± 0.05; 



1.01 ±0.01. 



(12) 
(13) 
(14) 

(15) 
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Table 1. Correlation function parameters of model samples 



Model 


N 


Step 


Train 


Tmax 


A-max 


A 2 i 


A 32 




CL.vor.00 


2920 


115 


81 


148 


0.073 








SC.cor.25 


1684 


120 


64 


146 


0.40 


127 


104 


120 


CL.rod.20 


3240 


115 


82 


138 


0.185 


119 


109 


116 


SC.rod.15 


3003 


115 


82 


141 


0.150 


118 


109 


115 


SC.nct.OO 


3452 


128 


96 


139 


0.45 


150 


118 


129 


SC.net.20 


2719 


115 


74 


142 


0.175 


119 


101 


116 


SC.mix.15 


4021 


115 


82 


136 


0.022 


121 


115 


116 




Figure 5. The correlation functions found for the whole cube (left panel) and for double conical samples (right panel). In both panels 
these functions are given for the regular net model and for the random supcrcluster model. Error corridors for the correlation function 
for the random supcrcluster model are plotted with long-dashed lines; error corridors for the regular net model with short dashed curves. 



As we see the most accurate determination of P comes from 
the last equation. 



5 DISCUSSION 

In this Section we continue the analysis of the correlation 
function and the power spectrum for different cluster distri- 
butions. Our main goal is to find which geometric properties 
of the distribution of superclusters can be detected on the 
basis of the correlation function and the power spectrum 
alone. We investigate the influence of the shape of the vol- 
ume occupied with clusters. The available observed samples 
are generally not cubical in shape, and the shape of sam- 
ple can, in principle, influence statistical properties of the 
samples. 



5.1 The correlation function 

The correlation functions with their error corridors and the 
respective power spectra with their error corridors are plot- 
ted for several models in Figs 2-4. On small scales all 
correlation functions have a large maximum which reflects 
the fact that our models have clustering properties. On large 
scales the correlation functions are different. They have one 
of three principal forms: they become flat on large scales 



directly from the maximum on small scales (random super- 
cluster model); have a minimum on intermediate scales fol- 
lowed by one secondary maximum, and then smoothly be- 
come flat (Voronoi model); or obey an oscillatory behaviour 
with alternating secondary maxima and minima of decreas- 
ing amplitude (models with a built-in regular structure of 
the type of a rectangular lattice). 

Geometric interpretation of the correlation function of 
the first type is simple: here superclusters are distributed 
randomly and peaks from numbers of objects at various sep- 
arations cancel each other out. 

Superclusters separated by large voids produce a cor- 
relation function with a minimum, whose location corre- 
sponds to the mean separation between superclusters and 
voids. The minimum is followed by a secondary maximum 
corresponding to the distance between superclusters across 
voids. On still larger scales there are no regularities in the 
supercluster-void network (voids are spaced randomly) , and 
the correlation function on very large scales is close to zero 
since superclusters at various separations cancel each other 
out as in the previous case. This distribution is realized in 
the Voronoi model. 

The interpretation of the oscillating correlation function 
is also straightforward. In all models which produce an os- 
cillating correlation function the distribution of clusters on 
large scales is quasi-regular in the sense that high-density 
regions form a fairly regular network with an approximately 
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constant grid length. In order to get regular oscillations the 
geometric structure must not be too regular. An example 
of very regular structure with superclusters located only in 
corners of the grid shows that the correlation function has 
a number of small maxima and minima without any regu- 
lar oscillations. Quasi-regular oscillations are generated in 
all models with certain regularity and certain irregularity: 
the supercluster-void network must have a constant overall 
scale but individual superclusters must be located flexibly 
within this network. The scale of the model is determined 
by the period of oscillations. To check the correspondence 
of the model to reality one can use small deviations of the 
correlation function from strict regularity expressed through 
parameters A21 and A32. As we see (Paper II), deviations 
of the observed correlation function from strict regularity 
are what is expected for models with constant scale and a 
certain scatter in individual superclusters positions. This co- 
incidence cannot be accidental, since a visual inspection of 
the distribution of real superclusters also shows the presence 
of such a network of superclusters and voids (Paper I). 

The basic difference between the regular and Voronoi 
models lies in the correlation of positions of high-density 
regions on large scales: in the regular model positions are 
correlated and in the Voronoi model not. 



5.2 Influence of the shape of the sample 

Real cluster samples are limited to a certain maximal dis- 
tance, and clusters in the galactic equatorial zone are not 
seen due to obscuration. Thus real observational samples 
generally have a double conical shape. To investigate the 
influence of such a sample shape we selected from cubical 
cluster samples double conical subsamples. These subsam- 
ples have a conical height equal to half the size of the total 
sample, and the equatorial belt (galactic latitudes less than 
±30°) is excluded. The volume of these subsamples is about 
one-quarter the volume of the whole cube, and they con- 
tain about one-quarter the number of clusters than is in the 
original sample. 

In Figure 5 we compare correlation functions of two 
models - the regular net model and the random superclus- 
ter model - calculated for the whole cubical volume, and for 
a double conical volume. The first model is regular and the 
corresponding correlation function is oscillating; the second 
model is not. This Figure shows that the oscillating proper- 
ties of the sample are preserved in the smaller double conical 
sample. Only the error corridor is larger in the smaller sam- 
ples, as expected. This similarity of the oscillating properties 
is due to the fact that the small samples are large enough 
to contain in sufficient quantity the structural elements re- 
sponsible for oscillations. 

5.3 The power spectrum 

Spectra for models are plotted in Figures 2-4. Power spec- 
tra have three regions, corresponding to inhomogeneities 
with high, medium, and low spatial frequency (or small, 
medium, and long wavelengths) . This division is given by the 
maximum of the spectrum. Our calculations show that on 
small wavelengths (large wavenumber k) spectra are rather 
similar. This result is expected as perturbations on these 
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Figure 6. The integral distribution of void radii for selected 
models. Models arc plotted as follows: solid line - Voronoi 
model CL.vor.00; dash-dotted line - random supcrclustcr model 
SC.ran.20; dashed line - superclusters on rods model SC. rod. 20; 
dash-dot-dot line - superclusters in net SC. net. 20; dotted line - 
superclusters in corners model SC. cor. 25. 

wavelengths are determined by the internal distribution of 
clusters within superclusters, and superclusters in our mod- 
els are generated in a similar way. 

Calculations by Frisch et al. (1995) show that inho- 
mogeneities on very large wavelengths do not influence the 
distribution of clusters in models, instead superclusters are 
modulated, i.e. in some regions they are very massive and 
in others poor. Regularity of the structure is given by power 
on intermediate wavelength. 

Analytic calculations confirm this result that models 
with only slightly different power spectra in the medium 
wavelength region may have very different correlation func- 
tions. The mixed model with quasi-regular network of rich 
superclusters and voids has a spectrum which differs from 
the spectrum of the random supercluster model only by 
about 25 % near the maximum. Differences between power 
spectra of the Voronoi and regular rod models are also small 
whereas correlation functions are very different. Such small 
differences can result from different geometries of the super- 
cluster population. These examples show that the correla- 
tion function is much more sensitive to the presence of a 
small regularity in the distribution of clusters than it is to 
the spectrum itself. 

5.4 The distribution of voids 

One feature of the correlation function is not yet explained. 
Why does the Voronoi model show a minimum and maxi- 
mum in the correlation function and the random superclus- 
ter model does not? One possible reason could be a differ- 
ence in the distribution of voids. In the random superclus- 
ter model voids can be very small, which is not the case for 
the Voronoi model. To check this possibility we extracted 
voids from our model cluster samples and calculated the 
integrated distribution of void radii. Here we used proce- 
dures described by Einasto, Einasto & Gramann (1989) and 
Einasto et al. (1991). Results for the void radii are shown in 
Figure 6. 

Our calculations show that for most of our models the 
distribution of void radii is very similar. A factor which 
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influences the distribution of void radii is the presence of 
weak filaments between superclusters. If clusters in high- 
density regions are distributed similarly in different models 
then their correlation functions also are similar, as they are, 
for instance, for models SC. cor. 25 and SC. net. 20. Void di- 
ameters in these models are, however, completely different, 
as seen in Figure 6. In models with no filaments between su- 
perclusters (SC. cor. 25) voids are much larger - their mean 
diameter corresponds to the diagonal of the net - whereas in 
all models where rods between corners are also populated, 
void diameters are determined by distances between these 
rods, i.e. by the scale of the net. 

We conclude that voids are determined by properties 
of small filaments and the correlation function by the dis- 
tribution of clusters in high-density regions (see also Paper 
I)- 



Table 2. 


Cosmic variance parameter 


Model 


(AO 




b 


CL.ran.OO 


1000 


0.014 


0.64 


SC.ran.OO 


2890 


0.029 


1.58 


f^T ...... nn 


2920 


0.020 


1.08 


SC.cor.25 


1684 


0.040 


1.62 


CL.rod.20 


3240 


0.009 


0.53 


SC.rod.15 


3003 


0.032 


1.75 


SC. net. 00 


3452 


0.017 


1.00 


SC.nct.20 


2719 


0.029 


1.54 


SC.mix.15 


4021 


0.017 


1.08 


DPS. 6 


2388 


0.020 


0.98 


CDM.61 


2588 


0.036 


1.83 



and for the cosmic variance: 



6 ERROR ESTIMATES 

The errors in the two- and three-point correlations are dis- 
cussed by Mo, Jing & Borner (1992). Errors in the correla- 
tion function depend on high-order correlations, which are 
different for different geometrical distributions of clusters. 
The difficulty of the determination of errors lies in the fact 
that they depend on the distribution of clusters in the en- 
semble, and there is no unique comparison ensemble of sam- 
ples available. Thus parameters describing the errors of the 
correlation function must be determined for all different ge- 
ometrical distributions of clusters. 

In the calculation of the correlation function the number 
of particles in the Poisson sample is taken as very large. In 
this case we can ignore errors of the Poisson sample and 
calculate the error of the correlation function from the error 
of the number of pairs (DD{r)). Mo et al. give a formula to 
calculate the error of (DD(r)} through moments of the two- 
point and three-point correlation function. The first terms 
have the form: 

a 2 DD (r) = {DD(r)) + ^{DD(r)) 2 . (16) 

Here N is the total number of clusters in the sample, and 
b is a constant which depends on the high-order correlation 
functions. The first term of the Eq. (16) is the Poisson error, 
due to random errors in sampling of galaxies; the second 
term is the cosmic error or variance, due to variations in the 
distribution of clusters in different parts of the Universe. 

For our models we have calculated errors in the cor- 
relation function from the scatter of different realizations 
of a particular model, and the Eq. (16) was used only as 
a convenient parametrisation to determine the value of the 
parameter b. As we see below this formula is quite accurate 
and gives results in good agreement with direct estimates. 
To see better the dependence of errors on the number of clus- 
ters in the sample we find approximate formulae for random 
and cosmic errors of the correlation function on large scales. 
On these scales the correlation function £(r) ~ 0, and thus 
{RR{r))(n 2 /n 2 R ) w (DD(r)). We get for the random error 
in the correlation function: 



0~£c 



N 



(18) 



a (P (r) : 



(17) 



The number of object pairs at separation r is propor- 
tional to the number of clusters squared, {DD(r)) ~ TV 2 , 
thus <7{p ~ A -1 . Random errors depend on the bin length 
dr and on the separation r. To reduce random errors we have 
smoothed the correlation function with a Gaussian kernel of 
dispersion a, which can be varied within certain limits. 

The cosmic variance is independent of the bin size and 
scale. It depends only on the geometry of the distribution of 
systems of clusters in the sample. The constant b increases 
with the size of the irregularities in this distribution. 

The influence of random and cosmic errors is illustrated 
in Figure 7. In the left panel we compare the effect of random 
and cosmic errors for samples of random clusters. The scat- 
ter of individual values in the unsmoothed correlation func- 
tion (circles) shows the effect of random errors, the overall 
deviation of open and filled circles from the mean correla- 
tion function (solid curve) is a visual demonstration of the 
effect of the cosmic variance. 

The influence of the smoothing length is studied in more 
detail in the accompanying paper devoted to the analysis of 
real cluster samples (Section 4.4 in Paper II). The number 
of clusters in observed samples is in some cases only sev- 
eral hundreds, thus random errors are larger. For observed 
samples a smoothing length about 15 h^ 1 Mpc reduces the 
random noise in practice but does not influence parameters 
of the correlation function. Our models are designed for com- 
parison with observations, thus we use a smoothing length 
15 /i -1 Mpc in most model samples too. 

In the right panel of Figure 7 we compare the error cor- 
ridors found from the scatter of individual realizations of 
the random supercluster model and calculated from the Eq. 
(16) or (18) (both formula give practically identical results 
on large scales). We see that the widths of both error corri- 
dors are approximately the same. This shows that the Eq. 
(18) is a good approximation of the cosmic variance. 

To derive the possible cosmic variance we have used not 
only geometrical toy models but also several A-body simu- 
lations of structure formation and evolution. These simula- 
tions were performed with 128 particles and 256 grid cells 
using a particle-mesh code. The size of the computational 
cube was taken to be L — 700 hT 1 Mpc w 6Ao, where Ao 
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Figure 7. Random and cosmic errors of the correlation function of model samples. In left panel we plot by open and filled circles values 
of the unsmoothcd correlation function for two realizations of random samples of 1000 clusters. The solid curve is the mean smoothed 
correlation function for ten realizations using a Gaussian kernel with dispersion a = 15 h~ 1 Mpc. Dashed lines around the smoothed 
curve show the error corridor calculated from the scatter of ten realizations. In the right panel we show the mean smoothed correlation 
function for the random supercluster model; dotted lines show the error corridor calculated from the Eq. (16), and dashed curves are for 
the error corridor calculated from the scatter of individual correlation functions for ten realization. 



is the characteristic scale of the supercluster-void network. 
Clusters of galaxies were selected by a procedure similar 
to the "friends-of-friends" algorithm. The lower limit of the 
cluster masses in the simulations is determined by the num- 
ber of clusters in the sample. The cube was divided into 
8 subcubes of half the size of the original cube; perturba- 
tions within subcubes can be considered as independent of 
each other. The correlation function of clusters of galax- 
ies was calculated for all subcubes, and the cosmic variance 
was found from the scatter of results of individual subcubes. 
The power spectrum for the sample of clusters of galaxies 
was calculated as for all other cluster samples separately for 
all subcubes, and then the mean value and its scatter were 
found. 

In simulations we used a CDM spectrum with the 
standard transfer function (Bond & Efstathiou 1984). This 
model is COBE normalized in amplitude at long wave- 
lengths, and it has one free parameter, T = Qh, depend- 
ing on the total density of the model £1 and the Hubble 
constant h. The model is designated as "CDM. 61", where 
6 designates the size of the computational box in units of 
the scale of the supercluster-void network, and 1 is for the 
standard high-density model (fi = 1). 

To estimate the value of the parameter b from the cos- 
mic variance Eq. (16) we found the mean la rms deviations 
of individual values from different realizations of the corre- 
lation function for the scale range r/L — 0.25 to r/L — 0.85. 
(Here L is the size of the computational box.) In this scale 
range rms deviations are approximately constant. The pa- 
rameter 6 is given in Table 2. 

The parameter b depends primarily on the regularity of 
the distribution on large scales. If there are no superclusters 
present (random cluster and clusters in rods models), the 
parameter is less than 1. In the Voronoi model and regular 
net model without shifts the error parameter has a value b 
1. In models with superclusters which have more freedom in 
their position (say, to move along rods) the parameter has 



a larger value 6 w 1.8. The CDM. 61 model also has a value 
of the parameter b close to this. In this case, however, the 
amplitude of oscillations is smaller than for the sample of 
clusters in high-density regions. A net model (superclusters 
in corners with additional clusters on rods) with scatter 15 
h^ 1 Mpc around the mean position is probably closest to 
the real case. 

To determine from a sample of objects the presence of 
oscillations of the correlation function, the cosmic error must 
be smaller than the amplitude of oscillations. If we adopt 
b — 1.6 and assume that oscillations have an amplitude 0.28 
(Paper II), then the sample must have at least 30 clusters. If 
we have a mixed population with an expected amplitude of 
oscillation 0.05, then such oscillations can be detected with 
confidence if the sample has at least 1000 clusters. We see 
that to detect oscillating properties it is useful first to isolate 
the population with highest level of oscillating nature. 



7 PROPERTIES OF THE CORRELATION 
FUNCTION ON SMALL SCALES 

Our main emphasis was the study of properties of the corre- 
lation function on large scales. Our models contain, however, 
information also on the correlation function on small scales. 
Thus it is of interest to have a look on small scales, too. 

As it is well known, on small scales the correlation func- 
tion can be approximated by a power law. This law can be 
interpreted in terms of the fractal distribution of galaxies. 
The fractal dimension is related to the power index (Szalay 
and Schramm 1985). The fractal description is the princi- 
pal geometric interpretation of the correlation function on 
small scales, and describes the shape of systems of galaxies. 
On small scales spheroidal systems of galaxies dominate, 
whereas on larger scales the main structural elements are 
filaments of galaxies (Einasto 1992). 

It is well known that clusters of galaxies have much 
larger value of ro than galaxies (Bahcall & Soneira 1983, 
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Klypin & Kopylov 1983) . This result is usually interpreted as 
an indication of the stronger clustering of clusters of galax- 
ies. However, Einasto et al. (1994) and Jaaniste et al. (1997) 
showed that the mean minor diameter of superclusters is 20 
h^ 1 Mpc. We have calculated the correlation length sepa- 
rately for clusters located in rich and poor superclusters, 
both in geometrical models and CDM simulations. As in 
the real case, supercluster richness was determined by the 
clustering analysis. Our calculations show that the correla- 
tion length is very different for these populations, about 17 
h^ 1 Mpc for clusters in poor superclusters and about 45 
h^ 1 Mpc for clusters in rich superclusters. Similar results 
come from the analysis of real clusters in Paper II. 

Thus, a correct geometric interpretation of the correla- 
tion length of clusters of galaxies is to say that systems of 
clusters of galaxies (superclusters) are larger than systems of 
galaxies (clusters) since the parameter ro depends primarily 
on the size of the respective systems. 



8 CONCLUSIONS 

Our study shows that the correlation function has the fol- 
lowing basic properties. 

On large scales the correlation function of clusters of 
galaxies depends on the geometry of the distribution of su- 
perclusters. If superclusters are located in a quasiregular 
net (supercluster-void network) then the correlation func- 
tion has an oscillatory behaviour. The period of oscillations 
is equal to the length of the step size of the net. The am- 
plitude of oscillations depends on the scatter of distances 
between superclusters in the net and on the presence of a 
more smoothly distributed population of clusters. 

Parameters of oscillation and the scale of the net can 
be determined from observations if the number of clusters 
is sufficient to get an error corridor of the correlation func- 
tion which is smaller than the amplitude of oscillations. For 
a quasiregular population this minimal number of objects 
is about 30, and for a heterogeneous population it is about 
1000, depending on the fraction of the quasiregular compo- 
nent in the whole population of clusters. 

If superclusters are located in the corners of a random 
cellular void network (Voronoi model) then the correlation 
function has a minimum, and its location corresponds to the 
mean separation between superclusters and voids. There fol- 
lows a secondary maximum corresponding to the clustering 
of clusters on opposite sides of voids. 

If the distribution of superclusters is random then the 
correlation function is featureless and becomes flat after the 
initial highly positive value. 

Oscillating properties of the correlation function are re- 
lated to the shape of the power spectrum near the maxi- 
mum of the spectrum. Oscillations occur only in the case 
when the spectrum has a sharp maximum and sudden tran- 
sition of the spectral index from a positive value on large 
wavelengths and a negative one on short wavelengths. 

Thus we see that the correlation function is indeed a 
useful statistic of the geometry of large scale structures. Ob- 
servationally, we have employed the large scale correlation 
function of a sample of Abell clusters (E97, Paper II); we 
have exploited these results in a study of the correlation 



function of Las Campanas Redshift Survey galaxies (Tucker 
et al. 1995, 1997). 
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